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I. INTRODUCTION 

Recently two orbitally excited narrow B s mesons were reported to be observed by the 
CDF collaboration with masses m B , = 5829.4 ± 0.7MeV and m B * = 5839.6 ± 0.7MeV 

. The DO collaboration confirmed the B* 2 state with mass m#* 2 = 5839.6 ± l.l(stat.) ± 
0.7(syst.)MeV subsequently, but their searching for a B sl signal gave inconclusive results 
with the available data set [2|. According to heavy quark effective theory (HQET) jjjj], 
the heavy quark spin decouples from the light degrees of freedom in the heavy quark 
itlq oo limit. The heavy-light mesons can be categorized according to the total angular 
momentum of the light degree of freedom j q . Combining j q with the spin of the heavy 
quark yields a doublet of heavy-light meson states j = j q ± 1/2. For 5*- wave ground-states 
with j q — 1/2 (L — 0), two states of negative parity form the H doublet (0~,1 - ). For 
P-wave excited states with j q = 1/2 or j q = 3/2 (L = 1), four states of positive parity form 
the S doublet (0 + , 1 + ) and the T doublet (l + ,2 + ). The two observed excited B s states 
above are considered as members of the T doublet in the bs system [lj]. 

The four P-wave excited bs states are often referred to as B**. If kinematically allowed, 
they are expected to decay into B*K, BK or both. The S doublet of bs system, (B* , B' sl ), 
are supposed to be broad because they decay through an S'-wave transition while the T 
doublet, (.B s i, ^2)) are narrow because they decay through a P>-wave transition |4j. This 
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is the case of B s i and B* 2 because their masses lie above the thresholds of B*K and BK 
and they can decay dominantly through these strong channels as is expected. But we don't 
know the masses of B* Q and B' sl experimentally. In fact, they are speculated to lie below 
the thresholds of B*K and BK [5]. Therefore, the channels (B* ,B' sl ) (B*,B) + K 
will be kinematically blocked and the isospin violating strong decays and electromagnetic 
decays will be their dominant decay modes. If so, they should be very narrow as well. This 
is quite similar to the situation of similarly excited D s mesons jsj. 

Lots of work have been done to study the properties of these states and their decays. 
Their masses had been predicted theoretically in Refs. 

SB 

, IsJ through different methods. 
The strong decays were investigated in Refs. 0, using different models. Motivated by 
the discovery of 5 s i(5830) and -E>* 2 (5840), we investigate semileptonic decays of B** into 
the H (D s (1968), £>*(2112)) and S (L> sJ (2317), £> sJ (2460)) doublets of cs mesons using 
QCD sum rules in the framework of HQET in this work. 

The QCD sum rule method 10], in combination with HQET, has been widely used 
to explore the properties and decays of heavy-light mesons and has been proved to be a 
powerful theoretical apparatus for studying heavy hadrons. In the past twenty years, study 
of semileptonic decay of the ground-state B and B s mesons by this method and quark 



models has been extensively explored 11 



-19j|. In our previous works, we have investigated 
the semileptonic decays of ground-state B mesons into some highly excited charmed meson 
doublets through this method 17]. Based on the general method of parametrizing the 
hadron matrix elements as proposed in Ref. [20J and the general interpolating currents for 
heavy-light mesons in HQET 2l|, it can be easily extended to the case of the semileptonic 
decays of excited B and B s mesons. To date, no work about the semileptonic decays of 
B** has been done with the same method. 

The remainder of this paper is organized as follows. After an introduction, we derive 
the formulas of weak current matrix elements and decay rates in HQET in Sec. [Ill In Sec. 
IIHI we deduce the three-point sum rules for the relevant universal form factors. Finally, 
section [IV] is devoted to numerical results and discussions. 



II. DECAY MATRIX ELEMENTS AND DECAY WIDTHS 

The semileptonic decay rate of a B meson transition into a charmed meson is deter- 
mined by the corresponding matrix elements of the weak axial-vector and vector currents 
(yv = crf^b and = 07^756) between them. The parametrizations of these hadronic 
matrix elements are important because it is not easy to calculated them from first princi- 
ples. They are often parametrized in terms of form factors which describe the momentum 
distribution among composites of a hadron. In HQET, the number of form factors is re- 
duced dramatically and this brings about much convenience. Especially at the leading 
order of the heavy quark expansion, the transition between certain low-lying heavy quark 
doublets can be parametrized by only one Isgur-Wise form factor for each process. This 
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can be done in different methods. In Ref . 22] , they finished the parameterizations in CQM 



model. Here we make use of the trace formalism given in Ref. 20] . 



A. B sl (B* s2 ) -> D s (D* s )ei7 and B sl (B* 2 ) -> D s0 (D' sl )£v decays 

The heavy-light meson doublets can be expressed conveniently as effective operators. 
For the decays B sl (B* 2 ) ->■ D s (D* s )iV and B sl (B* 2 ) -»■ /^^(-D^)^, one needs to introduce 
six field operators to represent the six heavy-light mesons that appear in these processes. 
The fields P v and P*^ v that annihilate members of the initial T doublet with four-velocity 
v are in the representation 

T* = 1 ~^{P^1, ~ \j\p v lsM - ^(7^ - «")]}, (1) 

where f — v ■ 7. The fields P, P* Pq, and P[* that annihilate members of the final H and 
S doublets with four-velocity v are, in form, 

Bv = ^[FZf-P'1*\ (2) 

and 

^ = ^[^75 + Po]. (3) 

At the leading order of the heavy quark expansion, the hadronic matrix elements of the 
weak current between states annihilated by fields in H v / and T v are 

t&Th® = ay)Tr{v' a H v/ TT v a }, (4) 

where is the heavy quark field in HQET, H v > = To-ff^Tb, and £(?/) is a universal Isgur- 
Wise function of the product of velocities y(= v-v'). For hadronic elements between states 
annihilated by fields in S v > and T v , they are similarly written as 

f$Th® = ay)Tr{v' a S vl TT v a }. (5) 

Here we should notice that each side of Eq. (j3J) and Eq. (jSJ) is understood to be inserted 
between corresponding initial B s and final D s states. The hadronic matrix elements of 
B s i(B* 2 ) — > D s (D*)iU can be calculated straightforwardly from the trace formalism (jlj) 
and are given as follows: 

wm v-mB.^)) = _ n ■ w _ ^ _ (y _ 2)vV , + 3vVv , f 

y /m Bsl m Ds V 6 

+ i(y + l)v p v' T £^], (6) 
^e')\(V-Ar\B sl{v ^ fl x + ^ + /v) _ 3u w 

y/m Bsl mD% V 
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-(y + 1)K - ""V - i(y - l)e^ T (v T + v' T ) 

+ 2iv'N p v' T e P(7pr - iv'Pv p v' T e°^] , (7) 

{Ds{v ' m T ~ A ^ B M ML =^y)v^ a ,[(y + iy - t/V - ^p<^1, (8) 

y/rn B . 2 m D . 

-iv' a e^ T (v T + v' T )}. (9) 

For the decays B s i(B* 2 ) — > D s q(D' s1 )£v, the corresponding hadronic matrix elements are 
calculated from Eq. ([5]) as follows: 

V m Bsi m D s0 V 6 

- t(y - l)tv«^*n (10) 

^m Bsl m D > Bi V 6 

+ (2/ - l)(v" + + % + 1)£ /3<7MT K - <) 

+ 2iv'»v p v' T e PapT - iv^v p v' T e appT }, (11) 

--((y)v'<*e a0 [(l - y)ge» + ^ + iv p v' T e^\, (12) 



(L> s0 K)l(^-^)1^ 2 (^e)) 



--((y)v' a e':e a p[g^v CT - g^v* - - 
-z^> r -<)]. (13) 



In these matrix elements, e a (e' a ) is the polarization vector of the initial (final) vector meson 
and e Q( g is the polarization tensor of the initial tensor meson, v is the velocity of the initial 
meson and v' is the velocity of the final meson in each process. (V — A)^ = C7 M (1 — 75)6 
is the weak current. Using the matrix elements fl6|)-f lT3|) . we derive the differential decay 
widths for these processes in terms of the Isgur-Wise functions as follows: 

d ^(B sl -+ D s tV) = &AVA ^ m S m\v ~ l)Hv + 1)*[(1 + rl)(2y - 1) 

-2 n (y 2 -y + l)], (14) 
dT G 2 F \V cb \ 2 m% m 3 n * 3 5 

-2r 2 (5y 2 + y + 2)}, (15) 

^(B sl -+ D s0 iV) = G2AVA l^ A \m\y - l) f (l/ + I) 1 [(1 + rl){2y + 1) 

-2r 3 (y 2 + ?/ + l)], (16) 
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rir G%\V c h\ 2 m%> m 3 n , . , 

d j-{B sl D' sl iu) LJ± ^ \C(y)\ 2 (y - i)% + i) f [(1 + rj)(7y - 1) 

-2r 4 (5y 2 -y + 2)], (17) 

a -{B* s2 D s £77) = C 3607r f 2 D ' K(y)l 2 (y - i) § (y + i) f [(i + *i)(*v + 1) 

-2r 5 (3y 2 + ?/ + l)], (18) 

^(s; 2 d s *^7) = 3607r 3 s2 ' lta)l 2 (y - i)% + W + r|)(iiy - 1) 

-2r 6 (7y 2 -y + 4)], (19) 

^(s; 2 d^) 3 6 J 3 S2 IC(y)l 2 (y - 1)% + i) f [(i + - 2 )(% - 1) 

-2r 7 (3y 2 -y + l)], (20) 
^(^; 2 -> ^) = ggjr fll K(?/)l% - l)*(v + l)*[(i + + i) 

-2r 8 (7y 2 + i/ + 4)], (21) 
where r,; (i — 1, 2, • • • , 8) is the ratio between the mass of the final D s meson and that 



of the initial B** meson in each process. For example, r\ = . The only unknown 
factors in these equations (jl4l) - (l2Tl) are £(y) and ((y), which need to be determined by 
nonperturbative methods. 



B. B sQ (B' sl ) -> D s (D* s )iu and B M {B' A ) ->• D s0 (D' sl )£v decays 

The hadronic matrix elements and differential decay widths of B sQ (B' sl ) — > D s (D*)iV 
and B s0 (B' sl ) — > D s0 (D' sl )£i> can be calculated by repeating the procedure in the previous 
subsection. The only change is that the initial doublet is now the S doublet in stead of 
the T doublet. So, at the leading order of the heavy quark expansion, the hadronic matrix 
elements of the weak current between states destroyed by fields in H v > and S v are 

h$Th® = X (y)Tr{H v ,TS v }, (22) 

while those between states annihilated by fields in S v i and S v are 

t$rh® = K(y)Tr{S v ,TS v }. (23) 

Using formula ff22l and (123]) . one can easily derive the hadronic matrix elements of the 
decays B s0 (B' sl ) -> D s (D*)£u and B s0 (B' sl ) -»■ D s0 (D' sl )iV as follows: 

MI'MIM ^y-/), (24) 
^m B *m Ds 
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(D* s (v> : e>)\(V-Ar\B* s0 (v) 



y/m B:o m D * 
(D s (v')\(V-Ar\B' sl (v,e) 



y /m B > al m Ds 
(D* s (v',e')\(V-Ar\B' sl (v,e) 



^/m B >m Dt 
(D s0 (v>)\(V - Ay\B* s0 (v] 



^m Blo m Da0 
(D' sl (v',e')\(V-Ar\B* s0 (v) 



^m Bl m Ki 
(D s0 (v',e')\(V-Ar\B' sl (v,e] 



y/m B > Bi m Ds0 
(D' sl (v',e')\(V-Ar\B' sl (v,e) 



y/m B ' ai m D ' a 



--X{yy;[-{V - W + v'^ + iv P v' T e^ T \, (25) 

--X(y)e a [(y - 1)<7 Q " - vN« - iv p v' T e a ^}, (26) 

=x(yKe a [g a °(v» - v'») + g^v' a - g^v" 

-ie^{v r -v' r )l (27) 

=«(y)(i/" + ^), (28) 

--K{y)e':[-{y + l)<f " + t/V + (29) 

--K{y)e a [-{y + 1)<T + t;V a + iv p v' T e a ^\, (30) 

= - K( y y;e a [g a °(v» + v'") - g^v' a - g a »v° 

+ ie^(v T + v' T )\, (31) 



The differential decay widths of these semileptonic processes are readily calculated from 
these matrix elements. They are given in terms of xiy) an d n(y) as follows: 



(B* s0 D s iU) 



dT 

dy 



G 2 F \V cb \ 2 m 2 B *m s Ds 



48tt 3 
G 2 F \V cb \ 2 m 2 B * 171%* 
48tt 3 

-2r 10 (4y 2 -y + l)], 



\ x (y)\ 2 (y-l)Hy + l)Hr 2 9 -iy, 



\ X (y)\ 2 (y-l)-2(y + l)2[(l + r 2 w )(5y-l) 



(B* s0 £> s0 ^) 



G 2 F \V ch \ 2 m 2 B * m 3 D 



sO U >>° 



dr /T3l 



48vr 3 

' 48^ 
- 2r l2 (Ay 2 + y + 1)], 

1447T 3 

-2r 13 (4y 2 -y+l)], 
G 2 F \V cb \ 2 m 2 B , m 3 D 



«(|/)| 2 (y-l)2(i/ + l)Mrii + l) 2 , 



K(y)| 2 (j/ - + 1)2 [(1 + r 2 12 )(5y + 1) 



| x ( y )| 2 ( y -l) 5 ( 2/ + l)2[(i + r 2 3 )( 5y -l) 



1447T 3 

- 2r u {8y 2 + y + 5)] 



| X (y)|%-l)2( y + 1)3[(1 +r 2 4 )(l 3y+ l) 



G 2 F \V ch \ 2 m 2 B , m 3 Ds0 



1447T 3 

2r 15 (4?/ 2 + ?/ + !)], 



«(y)| 2 (y- l)Hy + 1)^(1 + r 2 15 )(5y + l) 



(32) 

(33) 
(34) 

(35) 
(36) 
(37) 
(38) 
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^-(b' s1 D' sl ev) = 3l \<v)\ (y - + + ^)(i3y - 1) 

-2r 16 (% 2 -y + 5)]. (39) 

The definition of (z = 9, • • • , 16) is the same as rj (z = 1, • • • , 8) in the previous subsec- 
tion. The only unknown factors in these equations, (!32|) - (l39l) . now are x{y) an d n{y) which 
are also nonpurterbative quantities. In the following section, we will apply the QCD sum 
rule approach to estimate these Isgur-Wise functions £(y), C{y)> x(y) an d K (y)- 



III. FORM FACTORS FROM HQET SUM RULES 

A. and ((y) 

In order to apply QCD sum rules to study these heavy mesons, we must choose ap- 
propriate interpolating currents to represent these states. Here we adopt the interpolating 



currents proposed in Ref. [2l| based on the study of Bethe-Salpeter equations for heavy 



mesons in HQET. These currents JjPj/ for the states with the quantum numbers j, P, 
ji in HQET are proved to have nice properties and satisfy 

(Oi^K?W,^> = fpjAfS Pp/ S kfe r,^ (40) 

and 

z(0|T[j;^(x)j;f;;J'(0)]|0) = 8 jjf 8 PP ,8 Jtji (-lY Stf* ■ ••</;''' J dt5(x - vt)TL PJt (x), 

(41) 

in the limit vtlq — > oo, where ■r] ai '" otj is the polarization tensor for the state with spin j, 
= g a/3 — v a v 13 is the transverse metric tensor, and S means to symmetrize the indices 
and subtract the trace terms separately in the sets (ai • • • ay) and (/?i ••• /?,•). fpj t is a 
constant and Hpj t (x) is a function of x, both of which depend on P (the parity of the 
meson) and jt (the total angular momentum of the light part). 

Following the remarks given in Ref. [l^] , we take the interpolating currents that create 
the initial excited B s doublet (B s i, B* 2 ) as 

^3- 1 
,3/2 - v-vy -jhvlfsiD? - -' 



4,+,3/2 = H)\/ ih vl5 (D? - -7fpt)«, (42) 



4^3/2 = [ -^T a ^h vlt ,D tu s, (43) 

where Df = D a — v a (v ■ D) is the transverse component of the covariant derivative with 
respect to the velocity of the meson. The tensor T al3,fJ,1/ is used to symmetrize the indices 
and is given by 

T*** = \{gTi v + 9t l - \gfgr- (44) 
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For the final charm-strange mesons, the currents correspond to the ground states P s (1968) 
and P*(2112) with spin-parity (0~, 1~) are 

4,-,i/2 = hvlbS, (45) 
4-.1/2 = ^ (46) 



and the currents that create the (D sJ (2307), P s j(2460)) doublet with spin-parity (0 + , 1 

1 

7! 



are 

1 „ 



J o,+,i/2 = -^h v (-i)p t s, (47) 



With these currents, we can now estimate the Isgur-Wise functions £(y) and ((y) using 



QCD sum rules. We follow the same procedure as Ref. [12j to study the analytical 
properties of the three-point correlators: 



* 2 / tfxd*ze« k - k - z \0\T[J Q ^ 1/2 (x)J^ v \0)J« +fl/2 (z)\0) = S^w, u , y)Z%, (49) 
d^e 1 **'*"*^ = ^(w, (50) 



where ■* = h(v') , ~f ll h(v) and J^' " ^ = hty^j^jshfo). The variables fc(= P — msif) 

and fe (= P' — m c v') denote residual "off-shell" momenta of the initial and final meson 
states, respectively. For heavy quarks in bound states they are typically of order Aq C d 
and remain finite in the heavy quark limit. Hi(u;,u; ,y) (i = 1,2) are analytic functions in 
the "off-shell" energies u = 2v ■ k and u' = 2v' ■ k' with discontinuities for positive values 
of these variables. They also depend on the velocity transfer y = v ■ v ', which is fixed in a 
physical region. Cy* A and C'y* A are Lorentz structures. 

Following the standard QCD sum rule procedure the calculations of Hj(u;, u) , y) (i = 1, 2) 
are straightforward. Take (149!) as an example. First, we saturate Eq. (1491) with physical 
intermediate states in HQET and find the hadronic representation of the correlator as 
follows: 

, • v f-^f+^y) 

^i,hadron{u, w , y) = -—= — — t ; — + higher resonances, (51) 

(2A + 3 — u — ie){2A_ i — u —le) 

where /_ i and / + 3 are the decay constants defined in Eq. (l40j) . Apj l is the bounding en- 
ergy of the heavy meson with total parity P and angular momentum of the light part ji. 
Second, the function Hi (a;, a; ,y) can be approximated by a perturbative calculation sup- 
plemented by nonperturbative power corrections proportional to the vacuum condensates 
which are treated as phenomenological parameters. The perturbative contribution can be 
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represented by a double dispersion integral in v and v plus possible subtraction terms. So 
the theoretical expression for the correlator has the form 



/peri/ • \ 
dvdv Pl , ^ 7 ' - ; — + subtractions + S'[ on ' i (cu, oj',y), (52) 
[y — oj — %e)[y — oj — le) 

where S[ 0Tld (u;, oj , y) is the nonperturbative contribution containing vacuum condensates. 
The perturbative spectral density and the coefficients of vacuum condensations can be 
calculated straightforwardly. The other correlator fl50l) can be dealt with in the same way. 

Assuming quark-hadron duality, the contributions from higher resonances are usually 
approximated by the perturbative continuum above a threshold. Equating the phenomenon 
and theoretical representations, the contributions of higher resonances in the phe nomenon 
representation (IBTj) can be eliminated. Following the arguments in Refs. [12J, |23|, the 
perturbative and the hadronic spectral densities cannot be locally dual to each other; the 
necessary way to restore duality is to integrate the spectral densities over the "off- diagonal" 
variable V- = v — v , keeping the "diagonal" variable v+ = ^4^- fixed. It is in v + that 
the quark-hadron duality is assumed for the integrated spectral densities. The integration 
region can be expressed in terms of the variables z/_ and u + and we choose the triangular 



region defined by the bounds: < v + < oj c , — 2J^^v + < z/_ < 2y^q-j-z/ + . A double 
Borel transformation in u and oj is performed on both sides of the sum rules, in which for 



simplicity we take the Borel parameters equal 12j, |l5|, |l6j: T\ = T2 = 2T. It eliminates 
the substraction terms in the dispersion integral ( |52l) and improves the convergence of the 
operator product expansion (OPE) series. Our calculation is confined at the leading order 
of perturbation. Among the operators in the OPE series, only those with dimension D < 5 
are included. For the condensates of higher dimension (D > 5), their values are negligibly 
small and their contributions are suppressed by the double Borel transformation. So they 
can be safely omitted. Finally, we obtain the sum rules for the form factors £(y) and ((y) 
as follows: 



UJcl 



i{y)f-,i/2f + , V 2Z~ {K -< 1,2+K+ ' 3/2)IT =^r—T^ / dv + e--[Avl + 3m s (y + V)v\ 

87T [y + L) J 2ms 

R K 1 (ss> 2 (9sS<r ■ Gs) m s 
- 6m s {y + l)u+\ - -^rm s — [1 - — 

4y + 5 mg 1 y + 5 ,a, . , > 

+ "I! - T^ J " 3x2^ (y + l)2 ( V GG> ' (53) 



1 1 "4- 

C(y)/+,i/ 2 / + ,3/ 2 e- (A +^+ A +.3/ 2 )/^ = i6?r / du + e-^ K - 6m s (y + l)v. 



3 



167r 2 (y 

+ 27m 2 s (y + l)u 



2m a 

■i,.. 1 ,,,2i {g s scr-Gs) 1 m as 
1 4 l 3 2T' 

s 



10 



B. x(y) and n(y) 

The derivation of sum rules for x{y) and K {y) is J us t the same. The initial states are 
now members of the doublet (B* Q , B' sl ) with spin-parity (0 + ,l + ), so the interpolating 
currents one should use here are (1471) and ( l48~l) . The final states are the same as those 
in the previous subsection. Different from the previous subsection, the correlators whose 
analytical properties we study now are as follows: 



z 2 / ^^^^'--^^(OlTfJo^.^^J^'^O)^^,/^^^) = ~ 3 (cu,cu',y)£^ A , (55) 
i 2 j dWze^ fc '--^(0|T[J 0j+ , 1/2 (^)jg' u ' ) (0)J t )+jl/2 (z)|0) = E 4 (cu,u',y)C% A , (56) 



By repeating the procedure in subsection II II A\ we reach the sum rules for x(y) an d ^(y) 
in HQET as below: 



W C 3 



x(y)f-V2f + V2e- {A -^ +R +^ =(-^5)7^2 / du + e-^[u 3 + - 3m s (y + \)v\ 



1 2m 

+ Qm 2 (y + l)v+) + ^P-m s [3 + m 



T 



(g s sa-Gs) y + 7 

H — 2w — 1 — m s — — — 

24T 4T 

7 ?/ — 1 .a 



3 X 2 6 ?/ + r 7T 

CJ C 4 



GG), (57) 



+ 3m s y(y + - 3m 2 (y + l)(2y + 1] 

_(M-^) (4y _ 1)(1 _^ 

16 (y + l) 2 V 7T 1 K ' 

IV. NUMERICAL RESULTS AND DISCUSSIONS 

We now turn to the numerical evaluation of these sum rules and their phenomeno- 
logical implications. The input parameters are as follows. The mass of the initial and 
final heavy mesons are M Bsl = 5829.4MeV, M B > sl = 5765MeV, M Ds = 1968. 5MeV, 



M D * = 2112.3MeV, M Ds „ = 2317.8MeV, M D > = 2459.6MeV W, Mb* = 5839.7MeV, 
and Mb s0 = 5718MeV p]. For the QCD parameters entering the theoretical expres- 
sions, we take the standard values: (ss) = -0.8 x (0.24) 3 GeV 3 , (a s GG) = 0.04GeV 4 and 
ml = 0.8GeV 2 . The mass of the strange quark is m s = 150MeV. In addition, V c b = 0.04 
and Gf = 1.166 x 10~ 5 GeV~ 2 . The cutoff parameter is chosen as fj, — lGeV [25]. 

Let's evaluate the sum rule for £(y) numerically first. As we can see in Eq.( l53j) . two 
decay constants (.f+,3/2 and /- 1/2) and two bounding energies (A + 3 / 2 and A_ jl / 2 ) appear 
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on the left-hand side as input parameters. In order to obtain information of Isgur-Wise 
function £ (y) with less systematic uncertainty, we can remove these parameters by dividing 
the three-point sum rules with the square roots of relevant two-point sum rules, as many 
authors did 
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15 



16 



This can not only reduce the number of input parameters but also 
improve stabilities of the three-point sum rules. The two-point QCD sum rules we need 
here are 



18 



ft 



fl 



,3/2 



e 2A + , 3/2 /T 



64tt 2 



due t (z/ 4 + 2m s u 3 - Qm 2 u 2 - 12m\v) 
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(59) 



(60) 



After the division has been done, we obtain an expression for the £ (y) as a function of the 
Borel parameter T and the continuum thresholds U&, and Co>2- Imposing usual criteria 
for the upper and lower bounds of the Borel parameter, we found they have a common 
sum rule "window": 0.7GeV < T < 1.5GeV, which overlaps with those of the two-point 
sum rules ( 159]) and fl60|) . Notice that the Borel parameter in the sum rule for three-point 
correlator is twice the Borel parameter in the sum rules for the two-point correlators. In 
the evaluation we have taken 1.9GeV < lu < 2.4GeV 0, [3 and 2.8GeV < uo 2 < 3.2GeV 



251 ] . The regions of these continuum thresholds are fixed by analyzing the corresponding 



two-point sum rules. Following discussions in Refs. 12|, |23( , the upper limit uj c \ for u + in 



the region |[(y + 1) — \Jy 2 — l]wo ^ lo c i ^ |(^o + uj-i) is reasonable. So we can fix u c in 
the region 2.4GeV < u c \ < 2.8GeV. The results are showed in Fig. 1, in which we fix 
u = 2.2GeV and u 2 = 3.0GeV. The resulting curve for £(?/) can be well fitted by the 
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linear approximation 

£(</) = £(1)[1 - p\(y - 1)], £(1) = 0.81 ± 0.09, p\ = 0.83 ± 0.06. (61) 

The errors reflect the uncertainty due to u c and T. 

The numerical estimation of ((y) can be done in the same way. The two point-sum 
rules we need to remove low energy parameters /+,3/2, f+,1/2, ^+,3/2 and A +j i/ 2 from Eq. 



(PI) are Eq. ([59]) and [25 1 



/+,i/2 e_2A+ ' 1/2/T / dve-$(v A + 2m,z/ 3 - 6m 2 s u 2 - 12m 3 s u) 

1 o, >, 4m?, 3 , , m« ,a 



-m^>(l ~ f + l ji) + ^(ss) - ^GG). (62) 

The Isgur-Wise function ((y) finally appears as a function of the Borel parameter T and 
the continuum thresholds u C 2, uj\ and L/J%. The region of 0J2 is given above while ui\ is fixed 
in the interval 2.7GeV < u\ < 3.1GeV [25|]. The upper limit u C 2 for v + in the region 
\[{y + 1) — \fy 2 — l]wi ^ uj c2 ^ ^(cji + 0J2) is reasonable. So it can be fixed in the region 
2.7GeV < u C 2 < 3.1GeV. The results are shown in Fig. 2 where we fix u>i = 2.9GeV 
and U2 = 3.0GeV. As has been done to £{y), the curve of ((y) is also fitted by a linear 
approximation 

C(y) = C(l)[l - p\{y - 1)], C(l) = 0.085 ± 0.010, pi = 1.76 ± 0.05. (63) 

Using the linear approximates for £(y) and C(?/)> we can calculate the semileptonic decay 
rates of processes B s i(B* 2 ) —> D S (D*)£V and B sl (B* 2 ) — > D s0 (D' sl )iU. The maximal values 
of y for these semileptonic processes are given in Table [H By integrating the differential 
decay rates over the kinematic region 1 < y < y max , we get the decay widths of these 
semileptonic decay modes. Although the widths of B s \ and B* 2 have not yet been measured 



experimentally, they were estimated early in Ref. [26[ to be around 1 MeV. Theoretically, 
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TABLE I: The maximal value of y for each process: y ma x = (1 + rf)j1vi (i = 1, 2, • • • , 8). 





D s ei7 


D*JV 


D s0 £V 


D' sl £v 


B s i 


1.64951 


1.56105 


1.45633 


1.396 


B* s2 


1.65183 


1.56317 


1.4582 


1.39772 



their strong decays were investigated in Ref. [9(. As we know, the main decay modes of 
these excited B s mesons are strong decays. Therefore we can approximately take the strong 
decay widths as the total widths for an estimation of order of the branching ratios of these 
processes. In fact, the two-body strong decay widths of B s \ and B* 2 are computed to be 
98keV and 5MeV in Ref. [sj. Using these widths, we estimate the order of the branching 
ratios of the semileptonic decays (see Table ILT|) . Results of constituent quark meson (CQM) 
model in Ref. [22I are also shown there. 

TABLE II: Predictions for the decay widths and branching ratios of B s i(B* 2 ) — > D s (D*)£i7 and 
B sl (B* s2 ) ^ D s0 (D' sl )lV 



Decay mode Decay width T (GeV) BR one Result of Ref. [22] BR of Ref. [22] 



B s i ~ 


-> D s £v 


(2.6 ±0.6) 


x 10" 


-15 


~ 10" 


ii 


2.1 x 10" 


-15 


~ 10" 


ii 


B s i - 


-> D*£v 


(6.7 ± 1.5) 


x 10" 


-15 


~ 10" 


-11 


4.9 x 10" 


-15 


~ 10" 


-11 


B*s2~ 


-> D s iv 


(2.7 ±0.6) 


x 10" 


-15 


~ 10" 


-12 


2.1 x 10" 


-15 


~ 10" 


-13 


D* 

B s2 ~ 


-> D*£V 


(6.8 ± 1.5) 


x 10" 


-15 


~ 10" 


-12 


5.0 x 10" 


-15 


~ 10" 


-13 


Bsi - 


> Dsolv 


(1.0 ±0.3) 


x 10" 


-18 


~ 10" 


-11 


8.7 x 10" 


-20 


~ 10" 


-16 


Bsi ~ 


>D> sl £v 


(1.9 ±0.5) 


x 10" 


-18 


~ 10" 


-14 


1.0 x 10" 


-19 


~ 10" 


-16 


B*s2~ 


> Dsotv 


(6.5 ± 1.5) 


x 10" 


-19 


~ 10" 


-16 


5.6 x 10" 


-20 


~ 10" 


-18 


B*2- 


>D' sl £V 


(2.2 ±0.5) 


x 10" 


-18 


~ lQ- 


-15 


1.2 x 10" 


-19 


~ lQ- 


-17 



As we can see from Table HIl the decay widths and branching ratios of B s i(B* 2 ) — > 
D s (D*)iV quite agree with the results of Ref. [22| while there are obvious deviations in 
B s i(B* 2 ) — > D s q{D' s1 )Cm. Our results are about one order higher than theirs in the latter 
cases. After a simple derivation, the formulas of differential decay widths in the present 
paper and Ref. [22J are exactly the same ((14)-(21) in Sec. HTland (19) from Ref. 22]). So 
the reason for the differences must be that the Isgur-Wise function ((y) estimated in the 
present paper and the corresponding ((ui) in Ref. 22] are different from each other. In fact, 
the Isgur-Wise function £(y) here roughly agrees with £(u) in Ref. [22j while ((y) is quite 
different from ((uj) there. As mentioned above, the Isgur-Wise functions are estimated 
through QCD sum rules in this paper. The QCD sum rule method is a model independent 
method with its basis on first principle. It has been widely used for studying properties 
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1 1.1 1.2 1.3 1.4 1.5 1.6 1.7 



T (GeV) 



FIG. 3: (a) Dependence of x(v) on Borel parameter T at y = 1. (b) Prediction for the Isgur-Wise 
functions x(y) atT = lGeV. 



and decays of heavy-light mesons and proved to be a reliable method. In Ref. [22J], the 
authors calculated the Isgur-Wise functions from CQM model which is an intermediate 
approach, not as rigorous and general as that of the effective meson Lagrangian, but allows 
for a smaller number of input parameters In addition to that, the nonperturbative 
effect of QCD is systematically taken into account in the QCD sum rule method while in 
CQM model, it is just treated as a suppression of large light quark momentum in terms of 
damping factors or cutoffs in the loop momentum integral 14]. To this end, it is reasonable 
to say that our results, which come from first principle calculations, are more reliable than 
model dependent calculations. Note that both calculations in this paper and Ref. [22] are 
confined at the leading order of heavy quark expansion, an account of corrections of higher 
order may improve the results. 

Following the same way in which we estimate £(?/) and ((y) we can also calculate the 
three-point sum rules fl57"|) and f )58|) for x(y) and n{y) numerically. The decay constants and 
bounding energies that appear in fl57|) and (1581 as input parameters are /-,i/2, f+,1/2, A—,1/2 
and A +1 /2- The two-point sum rules we need to remove them are (1601) and (1621 . The regions 
of the continuum thresholds are 2.3GeV < u c s < 2.7GeV and 2.7GeV < u C 4 < 3.1GeV for 
the three-point sum rules for x(u) an d K{y), respectively. The results are showed in Fig. 3 
and Fig. 4 (In these figures, we fix u = 2.2GeV, uj\ = 2.9GeV, and u 2 = 3.0GeV). Here 
we note that the Isgur-Wise function n(y) which describes the transition between heavy 
mesons with the same quantum numbers (B s (0 + ) — > D s (0 + )) is approximately normalized 
as n(y = 1) = 1, which is implied by the heavy quark symmetry [3]. 

The resulting curves for x{y) an d n(y) can be similarly parametrized by the linear 
approximation 

X (y) = " P\{y - 1)], = -0.16 ± 0.03, p\ = 0.38 ± 0.05; (64) 



k{ V ) = k(1)[1 - P \{y - 1)], k(1) = 0.99 ± 0.09, p\ = 0.7 ± 0.1. 



(65) 
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0.4 0.6 0.8 1 1.2 1.4 1.6 ' 1 1.1 1.2 1.3 1.4 1.5 



T (GeV) y 

FIG. 4: (a) Dependence of n(y) on Borel parameter T at y = 1. (b) Prediction for the Isgur-Wise 
functions K,(y) at T = lGeV. 

Here we would like to give some remarks on the uncertainties of the sum rules ( |53|) . ( |54|) . 
( )57|) . and (|58|) . For the lack of information about higher resonances and continual state, 
the choices of the thresholds lead to the dominant uncertainties in the predictions for the 
Isgur-Wise functions. The input nonperturbative parameters, namely the values of the 
vacuum condensates, will also lead to errors to the form factors. In our calculation, only 
the errors due to the choices of sum rule windows are considered. 

The maximal value of y for each process of semileptonic B s q and B' sl decays is showed 
in Table [TTTI Using these maximal values and the forms of linear approximations for x(y) 

TABLE III: The maximal value of y for each process: y ma x = (1 + r i)/^ r i {} = 9> 10, • • • , 16). 





D s iv 


D*JV 


D s0 £V 


D' sl £V 


B.so 
B' sl 


1.62451 
1.63504 


1.53821 
1.54783 


1.43617 
1.44466 


1.37746 
1.38526 



and n(y), one can compute the decay widths of B sQ (B' sl ) — > D s (D* s )tP and B s0 (B' sl ) — >• 
D s q(D' s1 )£i>. Considering that the main decay modes of B s o and B' sl are isospin violating 
decays and radiative decays, they are supposed to have widths of about 100 keV One 
can then roughly estimate the branching ratios of these decays. All the results are presented 
in Table ED 

We would like to address that B s0 and B' sl are still missing in experiments. Their 
masses were theoretically estimated through various methods which gave quite different 
values The branching ratios given in Table [IV] are calculated on the assumption 

that B s q and B' sl lie below the thresholds of B*K and BK. If this assumption is not true, 
B s0 and B' sl can decay through these modes and they both own a width of hundreds of 
MeV. Then the branching ratios of the semileptonic decays in Table HVl should be 3 orders of 
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TABLE IV: Predictions for the decay widths and branching ratios of B s q(B' s1 ) — > D s {D* s )iv and 
B s0 (B' sl ) -> D s0 (D' sl )£U 



Decay 


mode 


Decay width T (GeV) 


BR 




one Result of Ref.[22j 


BR of Ref.[22] 


B s o —> 


D S £U 


(1.3 ±0.5) x 1(T 16 


~ 10~ 


12 


2.5 x 10- 14 


10- 9 ~ lO" 10 


B s o — > 




(1.5 ±0.5) x IO" 16 


~ 10~ 


12 


2.5 x 10- 15 


10- 9 ~ 10- 10 


B' sl ~> 


D s eu 


(6.6 ±2.3) x 10~ 17 


~ io- 


12 


1.2 x 10- 14 


lO- 9 ~ 10- 10 


B'si^ 


D* s t» 


(2.1 ±0.7) x 10~ 16 


~ io- 


12 


3.8 x 10- 14 


lO- 9 ~ 10- 10 


B s o — > 


D s0 eu 


(1.2 ±0.3) x IO" 14 


~ io- 


10 


1.6 x 10- 15 


io- 10 ~ io- 11 


B s o — > 




(3.4 ± 0.6) x IO" 14 


~ io- 


10 


1.5 x 10- 15 


lO- 9 ~ 10- 10 


B' sl ^ 




(4.0 ± 0.7) x IO" 14 


~ io- 


10 


4.9 x 10- 16 


io- 10 ~ io- 11 


B' sl ^ 


D' sl eu 


(1.1 ±0.2) x 10~ 13 


~ 10" 


-9 


2.1 x 1Q- 15 


IO- 9 ~ 10- 10 



magnitude lower. The present precision of the experimental measurement of the branching 



ratio of the B s mesons has reached up to 10~ 7 ~ 10~ 8 241 ] . As we can see from Table HT1 and 
Table IIVt our numerical results indicate that the decay widths of these semileptonic modes 
are tiny and therefore exclude the possibility of finding them in experiments. The decays 
B s0 (B' sl ) — >■ D s o{D' sl )tu may be expected to be measured in the future LHCb experiment. 
If the experimental results agree with our prediction, it can be supportive to the ordinary 
as meson explanation for D s j(2317) and _D s j(2460). 

In summary, we have performed a study of the semileptonic decays of B* s * mesons into 
low lying charmed-strange mesons within the framework of HQET. Two states out of the 
four P-wave excited B s mesons, B s \ and B* 2 , have been recently measured by CDF and 
DO Collaborations while the other two, B s q and B' sl , are still missing. We employ QCD 
sum rules to estimate the leading-order universal form factors describing the transitions of 
these B** states into low lying D s mesons, including D s (1968), D*(2112), £> s j(2317), and 
D s j(2460). The predicted branching ratios of these processes are prohibitively tiny. It is 
worth noting that the Isgur-Wise function n(y) which parameterizes the hadronic matrix 
element of the weak current between the S doublet of bs system and the same doublet of cs 
system approximately satisfies the normalization condition, k(jj = 1) = 1, which is implied 
by the heavy quark flavor symmetry at the leading order of the heavy quark expansion. 
Although the branching ratios we predicted are very small, some of them are expected 
to be possibly observed in the forthcoming LHCb experiment. A measurement of these 
processes can provide some information on the structure of the -D s j(2317) and -D s j(2460) 
mesons. 
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